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An  analytical  formula  of  the  emittance  of  a  held  emitter  is  given.  In  contrast  to  thermal  and 
photoemission,  such  a  formula  contains  complexity  due  to  the  multidimensional  nature  of  the 
source.  A  formulation  of  emittance  is  given  for  one-  and  three-dimensional  (3D)  held  emitters.  The 
3D  formulation  makes  use  of  the  point  charge  model  of  a  unit  cell  emitter  coupled  with  a  trajectory 
analysis  to  follow  electrons  to  an  evaluation  plane  where  emittance  is  determined.  The  single  tip 
theory  is  extended  to  an  array  and  the  resulting  theory  predicts  the  emittance  of  a  Spindt-type  square 
array  of  emitters  0.2  cm  on  a  side  producing  2000  A/ cm^  is  23  mmmrad.  Theory  compares 
favorably  with  experimental  measurements  in  the  literature  from  ungated  and  gated  sources.  The 
impacts  of  several  complications  are  estimated:  the  effects  of  a  gate  for  modulating  the  emitter;  the 
inhuence  of  space  charge  within  the  unit  cell  on  the  beam;  and  constraints  imposed  by  modulation 
frequency,  emitter  dimensions,  and  rise/fah  time  requirements  for  turning  a  beam  on  and  off,  as 


determined  by  the  array’s  RLC  characterization.  ©  , 
[doi:  10. 1063/1.3267288] 

INTRODUCTION 

Emittance  is  a  hgure  of  merit  for  judging  the  quality  of 
an  electron  beam  and  is  increasingly  important  as  larger  cur¬ 
rents  are  conhned  to  smaller  areas.  It  refers  to  the  tendency 
of  a  beam  of  particles  to  spread  in  cross  section  as  the  beam 
propagates  and  is  given  by  the  product  of  the  beam  size  and 
its  divergence.  The  emittance  due  to  the  cathode  itself  cannot 
be  undone  with  subsequent  beam  manipulations  (“beam  op¬ 
tics”)  and  therefore  knowledge  of  the  cathode  emittance  for 
various  electron  source  technologies  is  important.  Estimates 
of  the  emittance  of  a  thermal  emitter  and  a  photoemitter  exist 
in  the  literature.  It  is  the  purpose  of  this  work  to  provide  an 
account  of  held  emission  sources  and  is  motivated  by  appli¬ 
cations  that  seek  high  brightness  electron  sources.  Such  ap¬ 
plications  include  but  are  not  limited  to  advanced  particle 
accelerators  and  free  electron  lasers  (EELs),^  millimeter- 
wave  vacuum  electronic  (VE)  amplihers  and  terahertz 
devices,^’^  electron  beam  lithography,'^’^  transmission  elec¬ 
tron  microscopes,^  and  others. 

In  particular,  for  EELs,  the  electron  beam  must  be  fo¬ 
cused  within  the  laser  beam  optical  mode  within  the  wiggler 
in  order  for  the  laser  threshold  to  be  reached.  The  mode  is 
characterized  by  a  waist  radius  w  and  a  divergence  0,  the 
product  of  which  is  proportional  to  the  wavelength  of  the 
device,  or  w0=X/  it.  The  waist-divergence  product,  or  “emit¬ 
tance,”  must  be  such  that  the  beam  resides  within  the  optical 
mode.  Past  EELs  were  operated  at  wavelengths  of  tens  of 
microns,  and  so  the  emittance  of  the  electron  source  could  be 
large.  In  contrast,  recent  EELs  seek  to  operate  at  a  micron 
scale  or  smaller  (for  an  x-ray  EEL,  much  smaller)  and  so  the 
emittance  must  be  smaller. 
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Eor  micro-  and  millimeter-wave  power  amplihers  that 
use  thermionic  cathodes,  requiring  the  sources  to  produce  a 
current  density  characteristic  of  that  in  the  beam  tunnel 
would  render  the  lifetime  of  the  dispenser  cathode  unaccept¬ 
ably  short^  and  therefore  the  beam  transverse  dimensions  are 
compressed  after  emission.  Compared  to  an  ideal  or  zero- 
emittance  beam  /(e  =  0),  the  current  density  /(e)  in  the  beam 
tunnel  is  smaller  by  a  factor  that  depends  on  beam  radius, 
magnetic  held  strength  B,  accelerating  voltage  frequency 
/  (inversely  related  to  R^),  harmonic  number  n,  and  emit¬ 
tance  8  via  1  -[/(8)//(0)] oc  {nfs! RijBYVij.  B  and  are  con¬ 
strained  by  technology  and  platform,  higher  frequency  appli¬ 
cations  are  desired,  and  higher  harmonics  may  be  sought. 
Cathode  dimensions  for  past  microwave  amplihers  were 
large,  but  as  the  “power  density”  performance  metric  (the 
product  of  the  power  and  frequency  squared)  has  climbed, 
the  dimensions  of  the  beam  tunnel  have  shrunk,  and  so  a 
reduction  in  emittance  is  needed. 

Here,  a  methodology  to  evaluate  the  emittance  of  held 
emission  sources  is  developed.  Eield  emitters  have  sharp 
apexes,  therefore  the  usual  one-dimensional  (ID)  methods 
used  by  thermal  and  photoemission  theory  require  an  ex¬ 
panded  account  to  treat  curved  emitters  in  order  to  follow 
electron  trajectories  to  an  evaluation  plane  in  which  e  is 
calculated.  A  point  charge  model  is  provided  to  treat  rotation- 
ally  symmetric  three-dimensional  (3D)  structures.  The  effect 
of  combining  single  emitters  into  arrays  on  overall  emittance 
is  shown.  Complications  such  as  gating  and  space  charge  are 
considered.  Einahy,  the  modulation  of  actual  held  emitter 
array  (EEA)  sources  is  constrained  by  layout,  tip  geometry, 
and  modulation  frequency,  and  the  consequences  on  how  the 
arrays  are  operated  to  produce  bunched  beams  is  described. 
The  reader  is  cautioned  that  conventions  in  the  literature  (and 

903-1  ©  2010  American  Institute  of  Physics 


Author  complimentary  copy.  Redistribution  subject  to  AlP  license  or  copyright,  see  http://jap.aip.org/jap/copyright.jsp 


Report  Documentation  Page 

Form  Approved 

0MB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  0MB  control  number. 

1.  REPORT  DATE 

APR  2009 

3.  DATES  COVERED 

00-00-2009  to  00-00-2009 

4.  TITLE  AND  SUBTITLE 

Emittance  of  a  field  emission  electron  source 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Research  Laboratory, Code  6843,  Electronic  Science  and 

Technology  Division,4555  Overlook  Avenue  SW,Washington,DC,20375 

8.  PEREORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIEICATION  OE:  17.  LIMITATION  OE 

ARSTRAUT 

18.  NUMBER  19a.  NAME  OE 

OE  PAGES  RESPONSIBLE  PERSON 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  Sume  US 

unclassified  unclassified  unclassified  Report  (SAR) 

14 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


014903-2  Jensen  ef  a/. 


J.  Appl.  Phys.  107,  014903  (2010) 


respected  here)  mean  some  symbols  take  on  more  than  one 
meaning  depending  on  context:  subscripts  are  therefore  used, 
and  meanings  are  defined  in  the  tables. 


THE  EVALUATION  OF  EMITTANCE  FROM  A  PLANAR 
SURFACE 

Methodology 

The  normalized  rms  emittance  for  a  particle  beam  propa¬ 
gating  in  the  z  direction  is  defined  by^ 

e«.rms= — -  {xk^f,  (1) 

me 

where  x  is  position  and  hk^  is  momentum,  in  keeping  with 
quantum  mechanical  methods  to  be  used  below  which  rely 
on  conjugate  variables.  Expectation  values  are  given  by 
{0)  =  / (9/(r,k)JrJk,  where  /  is  the  normalized  distribution 
of  the  emitted  electrons:  thus,  {x^)  and  are  moments  of 
the  distribution  with  analogous  equations  along  the  other 
axes  y  and  z.  It  is  often  assumed  that  {xk^^  is  negligible,  but 
this  must  be  shown,  and  in  the  case  of  a  single  emitter,  it 
cannot  be  neglected.  Further,  f^D{k^fYij^E{y)\  where  the  z 
direction  is  normal  to  the  surface,  hk^  is  the  momentum  into 
the  barrier,  D  is  the  transmission  probability,  and  /pp  is  the 
Fermi-Dirac  distribution,  parameters  upon  which  they  de¬ 
pend  are  in  Table  I  (and  the  notation  respects  earlier  work.^) 
For  energy  parabolic  in  momentum,  usage  of  E  rather  than  k 
in  D  and  /  is  convenient.  Evaluation  of  the  moments  is  per¬ 
formed  in  the  xy-plane:  for  thermal  and  photoemission,  that 
plane  is  the  surface  of  the  cathode,  but  for  field  emission,  the 
location  of  that  plane  and  the  distribution  on  it  must  be  es¬ 
tablished.  First,  however,  thermal  and  photoemission  emit¬ 
tance  are  contrasted  with  field  emission  emittance  in  ID. 


Thermal  emission 

For  uniform  thermal  emission,  the  emitted  distribution  is 
independent  of  position,  and  so  (x^)  =  p^/2.  The  transmission 
probability  is  a  step  function  in  energy,  making  the  distribu¬ 
tion  approximated  by  the  thermal  tail,  exp(-ySp£’).  Current 
density  is  proportional  to  the  first  moment  of  the  longitudinal 
velocity  hkjm,  and  the  Richardson-Laue-Dushman  equa¬ 
tion  results.  The  cross  term  is  zero  as  it  contains  integrations 
over  even  and  odd  functions  in  k^.  The  well-known  expres¬ 
sion  for  thermal  emission  results  and  is 


rms(thermal)  = 


Pc 


(2) 


Predictions  based  on  Eq.  (2)  are  smaller  than  measurements. 
For  high  power  FELs,  bunched  beams  are  created  by 
“squeezing”  the  electron  beam  from  thermionic  cathodes 
producing  low  current  density  subject  to  comparatively  low 
gradient  because  of  space  charge  and  lifetime  concerns:  the 
acceleration  of  the  beam  is  therefore  not  rapid.  The  conse¬ 
quences  are  that  emittance  increases  by  a  factor  of  10-20 
over  the  ID  theory  predictions. 


Photoemission 

Photocathodes  switched  on  and  off  by  short  laser  pulses 
generate  electron  bunches  of  high  brightness  and  short  bunch 
length.  Although  not  as  straightforward  as  for  thermionic 
emission,  the  emittance  of  a  photocathode  has  been  recently 
developed.  Briefly,  the  energy  implicit  in  the  transmission 
probability  is  augmented  by  the  photon  energy.  Quantum  ef¬ 
ficiency  (such  as  current  density)  is  proportional  to  the  first 
moment  of  the  longitudinal  velocity  and  the  Fowler- 
Dubridge  relation  QE^{hci>- (pY  results.^’^^  However,  for 
emittance,  the  transverse  momentum  in  vacuum  {k^y  i^ust 
also  be  increased  by  the  photon  energy  ho)  (Dowell  et  al^ 
have  shown  that  differences  in  the  photoemission  emittance 
between  Refs.  11,  20,  and  21  are  due  to  the  mistaken  neglect 
of  photon  energy  in  k^  in  the  former),  and  it  can  be  shown 

en,rms(photo)  ^  —  I  )  •  (3) 


Field  emission 


In  contrast  to  thermal  emission,  the  transmission  prob¬ 
ability  for  field  emission  is  not  a  step  function  in  energy. 
Field  emission  from  metals  is  due  to  the  quantum  mechanical 
tunneling  of  electrons  through  a  surface  barrier  from  the 
Fermi  level  jx.  The  Fowler-Nordheim  transmission  probabil¬ 
ity  for  the  image  charge  barrier^’^^  is 

-  In  £»fn(^z)  = 

-Ei2m^ny){E,-ix),  (4) 


where  E^={hky  12m  is  the  “forward”  energy.  The  forms  of 
i;(y)  and  r(y)  in  Table  I  are  due  to  Forbes.^^’^^  The  first  mo¬ 
ment  of  the  longitudinal  velocity  hkjm  results  in  the 
Fowler-Nordheim  equation.  In  cylindrical  coordinates  for 
which  (kJ)  =  (/:^)/2,  using  Eq.  (4),  and  invoking  the  zero- 
temperature  limit  of  /pD, 

^^2  ^  ^  So'dkpik^sf^^kll TTkpdkp 
t'odkpik^)  jf^pl  TTkpdkp 


2ki 


hF 


(5) 


where  the  Fermi  momentum  kp  and  the  strong  exponential 
nature  of  D{k^)  are  used.  Therefore,  the  ID  field  emission 
emittance  formula  is 


Sn, rms(field)  =  p, 


hF 


i2m(py^hiy)  -  1 


-1/2 


(6) 


ID  comparisons 

Ratios  of  Eqs.  (3)  and  (6)  with  the  thermal  emittance 
[Eq.  (2)]  eliminate  the  cathode  radius  parameter.  The  three 
equations  are  valid  under  different  conditions,  the  Fowler- 
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TABLE  L  Terms  and  values.  Definition  and  values  of  common  terms.  “Cu”  refers  to  copper,  “Mo”  to  molyb¬ 
denum,  and  “Ba-W”  to  a  barium  monolayer  on  tungsten. 


Symbol 

Description 

Dehnition/value 

Unit 

Fundamental  constants 

h 

Planck’s  constant/2  tt 

0.658  211  89 

eVfs 

c 

Speed  of  light 

299.792  46 

nm/fs 

m 

Electron  mass 

510  998.9 

eV/c^ 

O^fs 

Fine  structure  constant 

1/137.036 

q 

Unit  charge 

1.602  176  5X10-19 

C 

Ub 

Boltzmann’s  constant 

1/11604.505 

eV/K 

eo 

Permittivity  of  free  space 

8.854187  8X  10-11 

F/m 

Distribution  function  and  emission  probability 

E 

Energy 

h^FUm 

eV 

p. 

Chemical  potential  (Fermi 

7[Cu] 

eV 

level) 

5.9  [Mo] 

kp 

Fermi  momentum  term 

hkp=^2mp 

1/nm 

O 

Work  function 

4.5  [Cu] 

2.1  [Ba-W] 

4.55  [Mo] 

eV 

T 

Temperature 

K 

Pt 

Inverse  temperature 

llkpT 

1/eV 

F 

Field  term 

q  X  electric  held 

eV/nm 

Q 

Image  charge  parameter 

a;f,^c/4=0.359  991  1 

eV  nm 

0 

Emission  barrier  above  p 

eV 

/fd(£) 

Fermi-Dirac  distribution 

{l+exp[^7<£-/i)]}-' 

D(K) 

Transmission  probability 

Ratio  of  transmitted  to  incident 

current  density  for  given  k^ 

Thermal,  field,  and  photoemission  terms 

^RLD 

Richardson  constant 

120.173  49 

A/K^  cm^ 

•^rld(^) 

Richardson-Laue-Dushman  equation 

exp(-</>/A:ar) 

A/cm^ 

Elliptical  integral  function 

Kj) 

Elliptical  integral  function 

^W^i+iy^d-inW) 

J 

Elliptical  integral  argument 

y=i4QFI<^ 

Pc 

Cathode  radius 

cm 

(O 

Photon  frequency 

rad/s 

Field  emitter  terms 

Fa 

Background  (asymptotic)  held 

eV/nm 

Pn 

Field  enhancement  factor 

Eq.  (36) 

Ftip 

Apex  held 

eV/nm 

a, 

Radius  of  emitter  apex 

nm 

Approx,  radius  of  emitter 

Base;  length  scale  of  PCM 

pm 

a,  a^i 

Apex  radius  factors 

Eqs.  (10)  and  (36) 

(dimensionless) 

Ps 

Field  enhancement  factor  due 

Eq.  (38) 

q/nm 

to  gate 

No 

No.  of  actual  emitters  on 

No.  of  emitters  =  (2/Vo +1)^ 

pm 

one  side  of  an  array  quadrant 

MV) 

Fraction  of  tips  emitting 

Approximately  linear  in  V 

N 

Number  of  active  emitters  on 
one  side  of  an  array  quadrant 

{2N+lfi  =  {2N^+lfiA{V) 

1 

Unit  cell  length  scale 

Comparable  to  d^  (Table  IV) 

pm 

L 

Array  length  scale 

L=(2N+l)l 

mm 

Nordheim  transmission  probability  is  not  valid  for  low  field, 
and  field  emitters,  in  particular,  are  not  planar:  comparing  the 
ID  emittance  formulas  for  the  same  parameters  is  therefore 
not  illuminating.  Instead,  evaluate  thermal)  for  T 

=  1300  K,  £„  j.jns(photo)  for  F=0.1  eV/nm,  and 
for  F=1  eV/nm  (all  of  which  are  representative  values)  but 


otherwise  use  copper  parameters.  For  a  266  and  213  nm 
drive  laser,  £  (photo) /£  (thermal)  is  then  1.27  and  2.25,  re¬ 
spectively,  whereas  £  (field) /£  (thermal)  is  3.31.  The  ratios 
reflect  that  the  full  width  at  half  maximum  of  a  photo  or  field 
emission  distribution  is  larger  than  the  width  of  the  thermal 
tail  of  a  thermal  emission  distribution  (see,  for  example. 
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Fig.  27  of  Ref.  25)  for  metals. 

There  are  caveats.  Photoemission  and  field  emission  are 
high  brightness  sources,  capable  of  in  excess  of  100  A/cm^ 
from  small  emission  areas,  and  therefore,  much  smaller  ac¬ 
tive  cathode  areas  are  needed.  Moreover,  improvements  in 
photoemission  emittance  are  more  a  consequence  of  superior 
engineering  in  forming  the  beam  than  advances  in  physics, 
e.g.,  rapidly  accelerating  the  beam  from  the  cathode  reduces 
emittance  but  requires  emission  to  be  bunched  or  “gated.” 
Therefore,  it  is  somewhat  misleading  to  consider  same- sized 
thermal,  photo,  and  field  emission  sources  or  imply  that  the 
beams  they  produce  are  directly  comparable. 


THE  EVALUATION  OF  EMITTANCE  OF  AN  ARRAY  OF 
CURVED  EMITTERS 

Field  emission  current  densities  are  on  the  order  of 
10^  A/cm^,  but  emission  sites  have  nanometer-scale  tip  ra¬ 
dii,  and  as  a  result,  the  active  emission  area  is  approximately 
50  nm^.  The  sharp  curvature  launches  electrons  with  a  trans¬ 
verse  momentum  component  from  the  start.  Placing  emitter 
tips  in  an  array  increases  the  area  over  which  emission  oc¬ 
curs.  Therefore,  the  emittance  of  an  array  of  field  emitters  is 
different  than  Eq.  (6):  to  calculate  the  emittance,  the  position 
and  velocity  of  all  electrons  crossing  an  evaluation  plane  are 
required  (in  contrast,  for  thermal  and  photoemission,  the 
evaluation  plane  was  the  cathode  surface). 

The  expectation  values  are  therefore  evaluated  in  a  man¬ 
ner  dissimilar  from  the  evaluation  of  thermal  emission,  in 
which  current  density  in  the  evaluation  plane  was  indepen¬ 
dent  of  position.  Electron  trajectories  are  followed  from  the 
emitter  tip  to  the  evaluation  plane  and  are  calculated  using 
the  point  charge  model  (PCM).  By  such  means,  the  crossing 
point  will  be  related  to  the  starting  point  on  the  emitter  sur¬ 
face.  Scaling  terms  that  allow  different  sized  geometries  to 
be  compared  based  on  a  prolate-spheroidal  model  (PSM)  of 
an  emitter  can  be  found,  and  we  do  so,  but  caution  that  the 
PSM  is  not  a  substitute  for  a  model  such  as  the  PCM  to 
deduce  trajectories. 


Single  ungated  emitter  in  the  PSM 
Position  and  veiocity  in  evaiuation  piane 

If  inertial  effects  are  ignored  (a  massless  approximation), 
then  field  lines  overlap  trajectory  lines.  The  PSM^^~^^  pro¬ 
vides  field  lines  analytically.  Dimensionless  coordinate  sys¬ 
tems  as  in  Table  II  make  connections  with  the  PCM  trans¬ 
parent,  allowing  for  the  usage  of  scaling  arguments,  and  are 
therefore  used.  To  relate  to  physical  systems,  “dimension¬ 
less”  coordinates  can  be  expressed  in  the  “units”  of  Table  III, 
although  dimensioned  terms  shall  be  explicitly  given  when 
needed. 

The  prolate  spheroidal  (alternately,  “ellipsoidal”)  coordi¬ 
nate  system  is  related  to  the  cylindrical  coordinate  system 
{p,z)  by 

p  =  sinh  a  sin  fS, 


TABLE  IL  Coordinate  systems  and  relationships.  Relationships  between 
cylindrical  (specihes  trajectories),  prolate  spheroidal  (specifies  emitter  sur¬ 
face)  and  spherical  (specifies  spherical  emitter  apex)  coordinates. 


Symbol 

Description 

Definition/relation 

P,  z 

Cylindrical  coordinates 

a,  p 

Prolate-spheroidal  coordinates 

p=sinh  a  sin  p 
z= cosh  a  cos  p 

(Ps’Z,) 

Emitter  surface 

is  constant 

iPh,Zh) 

Evaluation  plane 

Zh  is  constant 

Zh 

Height  of  evaluation  plane 

cosh  a}^(p)=Zh^ cos  p 

Zo 

Height  of  emitter 

Zsi^) 

R 

Ratio  of  semimajor  to  semiminor  axes  for 
ellipsoid  representing  emitter 

coth 

V 

Ratio  of  evaluation  plane  to  emitter  height 

Zh!  Zq 

a,  6 

Radius  and  polar  angle  of  equivalent 

Ps^a  sin  6 

sphere  representing  ellipsoid  apex 

Zs-Zo-a{l-cos  0) 

for  ps<a 

^Zo-p^s^2 

z  =  cosh  a  cos  p. 


(7) 


where  a  and  p  are  the  angular  coordinates,  and  because  p 
and  z  are  dimensionless,  the  usual  dimensioned  coefficient 
on  the  right  hand  side  is  omitted.  The  emitter  surface  is  de¬ 
fined  by  a  constant  a=a^,  where  the  5'  subscript  designates 
“surface.”  Eield  lines  are  defined  by  constant  The  evalu¬ 
ation  plane  is  defined  by  z=Zh-  The  relationship  between  pro¬ 
late  spheroidal  surface  and  evaluation  plane  coordinates 
(along  side  a  polar  coordinate  {a ,  0)  equivalent  sphere  sur¬ 
face)  is  shown  in  Eig.  1.  Also  shown  is  the  field  line  (red  or 
curved)  compared  to  a  massive  particle  trajectory  (blue  or 
straight). 

The  field  lines  originate  on  the  surface  at  (p^,z^)  and 
cross  the  evaluation  plane  at  (p/^ ,  Zh)  •  Therefore,  the  relation¬ 
ship  between  p^  and  p/^  is 


Ph 


sinh  a/j 

- Ps^ 

sinh 


(8) 


where,  because  Zh  is  constant,  is  a  function  of  p  and 
therefore  p^,  and  because  is  a  constant,  z^  is  a  function  of 
(5  only.  If  the  ratio  of  the  semimajor  to  the  semiminor  axis  of 
the  ellipsoid  is  R,  and  the  ratio  of  height  of  the  evaluation 
plane  to  height  of  the  emitter  is  rj=Zhl Zs{^)  —  Zhl Zq,  then 
very  close  to  the  emitter  apex  where  emission  is  dominant, 
using  notation  summarized  in  Table  II,  we  find 


TABLE  III.  Dimensionless  units  conversion  factors.  “Dimensionless”  pa¬ 
rameters  are  understood  to  be  in  the  following  units  based  on  the  choice  of 
length  and  energy  scales.  The  “Representative  value”  column  expresses 
quantities  in  SI  units  if  flo=l  and  F^=10^  MV/m. 


Term  Definition/relation  Representative  value 


Length 

Energy 

Mass 

Velocity 

Time 

Wave  number 


f’»ao 


(mao/2F„)‘'" 
kfj  =  {ImF^ao) 


1  yUm 

10  eV 

5  110  998.9  eV/c^ 
1  875  537  m/s 
0.533  180  5  ps 
16.200  88  nm-i 
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FIG.  1.  (Color  online)  Geometry  of  an  ellipsoidal  field  emitter,  showing 
cylindrical  coordinates  along  the  emitter  surface  (ps,Zs)  and  along  the  evalu¬ 
ation  plane  {ph,Zh)-  The  equivalent  sphere  characterized  by  {a,0)  is  also 
shown.  The  red  (curved)  line  from  the  surface  indicates  a  field  line,  whereas 
the  blue  (straight)  line  shows  the  trajectory  of  a  massive  particle  emitted 
with  an  initial  velocity. 


{xk^)=^aQkQ{pkp),  (13) 

where  the  and  terms  of  Table  III  handle  units,  the  h 
subscript  on  p  is  suppressed,  and  evaluation  in  the  Zh  plane  is 
understood.  Therefore, 

fi 

^n,rms I  single  emitter  “  )(^p)  ~  (P^p)  )  •  (14) 

The  function  P{ph)  is  a  ratio  of  current  density  at  to 
current  density  on  axis.  Therefore,  there  are  two  consider¬ 
ations  to  P\  the  magnitude  of  the  current  density  at  the  emit¬ 
ter  surface  and  how  it  declines  through  beam  expansion  at 
the  evaluation  plane.  At  the  surface,  the  emission  current 
is  related  to  the  field  F(p^)  via  the  Fowler-Nordheim  equa¬ 
tion:  we  use  the  Murphy  and  Good  (MG)  form^’^^  with  the 
Forbes  approximation  to  i;(y)  and  the  log  terms  in 

u(y)  changes  the  MG  form  to 


Ph  “  Ps^izl- 1){R^ -  . 


(9) 


When  inertial  effects  are  included,  p^l Ps  will  differ  from  Eq. 
(9).  From  Eq.  (8)  and  the  definition  of  apex  radius. 


«  =  [-  ^lp=o]  ‘ 


cosh(a^) 
sinh^(a^)  ’ 


R{a)  = 


1/2 


(10) 


The  term  in  the  radical  of  Eq.  (9)  dominates  for  Spindt- 
type  emitters  so  that  R^-\  goes  as  1/(2  when  a<\. 

The  initial  velocity  is  approximated  by  the  relation¬ 
ship  between  the  current  density  and  the  number  density 
Pg  using  Je=qPe^e  is  directed  normal  to  the  emitting 
surface.  To  leading  order,  the  magnitude  of  the  velocity  is 


v^  =  ■ 


QPe 


fik 

5^D{k,)f^^{\Pdk 

m 

jD(k,)f^j,{k)dk 


hkp 


m 


(11) 


as  a  consequence  of  the  strongly  peaked  nature  of  the  inte¬ 
grand  in  k^.  The  final  velocity  magnitude  is  obtained  from  the 
initial  velocity  and  the  potential  difference  between  the  loca¬ 
tion  on  the  evaluation  plane  and  the  emitter  surface. 


Distribution  of  eiectrons  in  the  evaiuation  piane 


The  symmetry  of  the  unit  cell  emitter  implies  that  the 
distribution  of  electrons  will  depend  only  on  the  cylindrical 
coordinate  p  and  not  on  the  angle  cp.  Denote  the  distribution 
of  electrons  in  the  evaluation  plane  by  P{p}).  Eor  a  single 
emitter,  x=aQp  cos  cp  and  so 


,  2n,_  2/ lo^dcpjQjp  cos  <pfp{p)pdp\ 

Sl-dq>SlP{p)pdp  I 


(12) 


and  similarly  with  kj^=kQkpCOS  cp  it  follows 


J^{F)  =  AF^-'^  exp(-  B/F) ,  (15) 

where  A,  B,  and  k  are  given  by 


(16) 

and  where  B  is  equivalent  to  in  Eq.  (40)  of  Ref.  9, 

6065,  and  t(yj  =  l  +  (l/6e)  =  1.061  [y^,  optimizes  a 
quadratic  fit  to  v{y)  for  metals^®].  For  copper,  k=0.772  81, 
B=65.207  eV/nm,  and  A  =  62.182  A/cm2(cm/eV)2-^ 

The  current  density  in  the  Zh  plane  is  related  to  the 
current  density  on  the  surface  by  conservation  of  particle 
number,  for  which 


K  = 


82  2m 
9/i  V  O  ’ 


A  = 


167r^^d>r(yJ^ 


B  =  — ^/2ma>^ 
U 


A2  =  2t7Pj 


2 

~  ^^Ph'^hiPh)^Ph^ 


(17) 


where  A  2  is  the  emitted  charge  in  unit  time  in  a  ribbon  of 
width  Ap^  on  the  surface  and  Ap^  in  the  Zh  plane.  Close  to 
the  emitter  apex,  the  surface  is  given  by  ZsiPs)  ^ Zq- P^ / 2a  so 
that  dZs! dps^-Ps! CL.  Solving  for  then  identifies  P{pi)  in 
Eq.  (12)  as 


PiPh)  = 


Jhiph) 

4(0) 


1 4(^’tip)  /■ 


(18) 


The  field  at  the  apex  is  E(p^=0)=4ip.  Along  the  spherical 
apex,  F  falls  as  cos  6  [Table  III  and  the  discussion  fol¬ 
lowing  Eq.  (36)  below]  so  that 
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JelFiPs)] 

JeiFtip) 


=  exp 


=  exp 


1 

cos  eip,) 


B 


2F,^a 


(19) 


Therefore,  P(ph)  in  Eq.  (12)  is  approximated  by  a  Gaussian 
in  the  evaluation  plane,  or 

PiPh)  ==  exp(-  Cpl), 


_ B _ B 

C  CpsM  -  \)R^  +  1]  ~  -  1)  ’ 

(20) 

where  the  tilde  indicates  that  although  C  is  similar  CpsM. 
latter  is  only  approximate  and  so  C  is  to  be  numerically 
found  from  the  PCM. 

In  the  evaluation  plane,  the  radial  velocity  changes  with 
distance  from  the  axis.  A  good  approximation  is 

kpiPh)^  yiph-y^pl-  (21) 

The  y  coefficients  are  determined  numerically,  although  a 
crude  approximation  for  yi  alone  may  be  made  by  assuming 
motion  in  a  constant  field  near  the  axis,  or 


1 

1/2 

Lip 

Fo^o 

.F^Zoiv-  1)_ 

(22) 

where,  as  with  Eq.  (20),  the  right  hand  side  is  only  approxi¬ 
mate.  Using  Eq.  (21)  it  is  found 

<p^)<^p>  =  <p^>[r?(p^>  +  2rir3(p'^)  +  7^(p®>]> 

ipkp)^  =  T^iipY  +  2yi  ysip^Kp'^}  +  ^^(pV-  (23) 

Making  use  of  the  relation 

f  p^"P(p)pdp=  (  p^"e~^P^pdp=  (24) 

Jo  Jo 


where  r(^)  is  the  gamma  function,  the  emittance  for  a  single 
emitter  is  then  found  to  be 


^n,rms  I  single  emitter 


^|2) 


^0^0 


me 


(25) 


field  is  good,  so  extrapolating  to  other  conditions  or  configu¬ 
rations  from  a  numerical  baseline  example  is  possible. 


Multiple  ungated  emitters 
Periodicity  of  an  array 

Assume  all  field  emitters  of  an  array  are  geometrically 
identical  and  operate  under  the  same  conditions — the  ap¬ 
proximations  are  not  strictly  true  because  emitters  vary 
among  themselves  and  those  at  the  boundary  of  the  emis¬ 
sion  area  will  experience  a  higher  field.  Cartesian  coordi¬ 
nates  are  required  for  a  rectangular  array.  The  Gaussian  na¬ 
ture  of  P(ph)  allows  for  a  simplification  in  the  evaluation  of 
quantities  in  one  Cartesian  dimension  v,  such  that  the  mean 
value  of  a  function  0{x)  is 


/  0(x)f(x,kjdxdk^ 


(26) 


where  f(x ,  k)  is  constructed  from  P(p)  for  the  spatial  distri¬ 
bution  of  emitters,  but  for  which  an  explicit  form  is  unnec¬ 
essary.  All  that  is  important  is  that  for  an  array,  /  will  be 
periodic:  if  the  unit  cell  has  dimensions  /,  then  f(x-\-jl) 
=f(x)  for  integer  j.  If  the  length  of  the  side  of  an  array  is  L 
then 


rLI2 

J-L/2 


L/2  ^  rl/2 

f(x)dx^  2  I  fix  +  jl)dx 

i=-N  J  -I 


j=-N  J  -1/2 


rl/2 

)J  . 

J  -1/2 


(2A-H1)  f{x)dx. 


(27) 


where  (2A-i-  1)/=L  and  (2A-i- 1)^  is  the  number  of  emitters  in 
a  square.  Thus, 


rL/2 

J-L/2 


L/2  ^  rl/2 

x^f{x)dx  =  2  I 

j=-N  J  -1/2 


(x  +  jlff(x  +  jl)dx 


rl/2 

:(2A-Fl)j  x^f(x)dx 

J -1/2 

4ifyrfix)dx. 

\  i=-N  /  J  -1/2 


(28) 


where  the  y^  terms  cancel.  Using  a  generic  y^-Tl  and 
aQkQ=^X  10^  from  the  trajectory  analysis  (below),  and  Cps^ 
indicates  that  £„  — 360^/mc.  Therefore,  single  emitters 

are  exceptionally  low  emittance  and  high  brightness  sources. 
Eor  arrays  of  emitters,  the  yi  terms  do  not  cancel,  the  cross 
term  is  negligible  if  the  number  of  emitters  is  large,  and  the 
array  emittance  is  much  larger. 

It  must  be  emphasized  that  the  PSM  is  only  an  approxi¬ 
mation  for  the  shape  of  a  Spindt-type  emitter  and,  more  im¬ 
portantly,  because  it  neglects  electron  inertial  effects,  it  is 
increasingly  inaccurate  in  estimating  p/^  and  kp{pf)  as  the 
evaluation  plane  is  further  removed  from  the  apex:  if  77 >  1 
then  PSM  predicts  that  kpipf^)  is  vanishing,  which  is  incor¬ 
rect.  However,  the  parametric  dependence  on  geometry  and 


The  summation  on  the  second  line  is  analytic.  Therefore, 
(x2)  =  |yv(yv+i)/2  +  (x\, 


!^!_l2X^f{x)dx 

i-yix)dx  ’ 


(29) 


where  the  second  line  refers  to  the  unit  cell  or  single  emitter 
and  is  so  indicated  by  the  o  subscript.  Eor  field  emitters,  / 
vanishes  at  the  boundaries  //2,  but  the  thermal  result  is  re¬ 
covered  if  f=  1  (uniform  emission),  for  which  Eq.  (29)  be¬ 
comes  {x^)  =  [4N{N-\-l)-\-l]P/3=L^/12.  Eor  arrays,  the  N 
terms  dominate  the  unit  cell  term.  Similarly,  will  exhibit 
the  same  periodicity  as  /,  and  so 
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rm  ^  rll2 

I  0{x)f{x)dx  =2  I  0{x+  jl)f{x  +  jl)dx 

J  -L/2  i=-N  J  -I 


j=-N  ^  -112 

rl/2 

=  {2N  +  1)  J  0(x)f{x)dx. 

J -112 


(30) 


Therefore,  {k^)o=k^{kPo  as  long  as  P{II2)  is  negligible.  An 
analogous  derivation  shows  that  {xk^  =  {xk^o' 


Array  emittance 


Ignoring  boundary  terms  reflects  that  the  array  is  infinite. 
A  systematic  study  by  Rhee  et  al?^  for  finite  sized  arrays 
obtains  a  similar  expression  to  ours  in  their  case  1  (parallel 
beam)  example:  their  geometric  factor  Gji  is  equivalent  to 
(/^/3)(1-hA“0(A-h1/2)^  and  is  a  factor  of  1-h(1/A)  larger 
than  ours.  For  N  sufficiently  large,  as  with  usual  arrays,  the 
difference  is  negligible.  A  finite  sized  array  experiences 
space  charge  forces  that  cause  emission  from  the  edge  of  the 
array  to  bend  away  from  the  center.  This  is  analogous  to  the 
focal  length  parameter  of  Rhee  being  finite.  An  infinite  focal 
length  limit,  as  presumed  by  the  periodicity  assumption, 
means  that  the  beamlets  are  parallel:  if  they  were  not,  addi¬ 
tional  terms  of  order  \l{2N+\)  and  \I{2N+\Y  arise.  Using 
(x^)~A(A-i-l)/^/3  and  {xk^=kQaQ{pk^ !  2,  the  emittance  for 
an  array  is 

^n,rms  I  array 

-  —  -  Uoao)"<PV^  (31) 

me  y  6  ^4 

Observe  that  {pk p)^- {yi / CY  is  small  compared  to 

SO  that  the  cross  term  is  negligible,  and  so  to  leading 

order, 

1/2 

(32) 

Using  Eq.  (32)  requires  attention  to  three  subtleties:  first, 
and  C  are  evaluated  numerically  although  PSM  can  give  or¬ 
der  of  magnitude  estimates;  second,  I  is  a  length  scale  char¬ 
acteristic  of  the  unit  cell  of  a  single  emitter,  but  is  not  iden¬ 
tical  to  the  tip-to-tip  separation  unless  all  emitters  are  active; 
third,  (2A-i-l)^  is  the  number  of  active  emitters. 


^n.rmslarray 


fl 

y\ — (w 

me 


N{N+  1) 
6C 


The  PCM  and  trajectory  parameters 
Mono-  and  dipole  potentials 

The  PCM,  introduced  in  a  previous  study,  analytically 
gives  field  enhancement  and  tip  radii  parameters.  A  line  of 
charges  is  placed  such  that  the  resulting  zero  equipotential 
line  approximates  the  shape  of  an  emitting  structure.  Previ¬ 
ously,  the  PCM  stacked  charges  of  one  sign  above  the  z=0 
plane  to  approximate  the  shape  of  the  emitter  created  by 
melting  and  field  formation.  The  potential  then  is  given  by 

y(p,z)-F„aoyj-,-  , 

\aQ  Gq/ 


FIG.  2.  (Color  online)  Comparison  of  mono  (blue  or  cusplike)  and  dipole 
(red  or  ellipsoidal-like)  surfaces  in  the  PCM  for  the  parameters  r=0.5  and 
n  =  6. 


Vnip,z)  =  -z+{p^  +  z^)  +  ■.2-11/2’  (33) 

j=l  IP  +  \Z-Zj)  J 

where  Gq  is  a  (dimensioned)  characteristic  length  of  the  base 
and  (p,2)  are  the  laboratory  coordinates.  The  Xj  are  deter¬ 
mined  by  the  boundary  condition  y„(0,z„+i)  =  0,  which  re¬ 
sults  in  a  matrix  equation  that  can  be  inverted  to  recover  the 
Xj' s.  Importantly,  p  and  z,  and  likewise  y„,  are  dimension¬ 
less,  allowing  for  scaling  arguments  to  be  made.  The  point 
charge  locations  are  defined  by 


n-l 

;=i 


1-r" 

1  -  r  ’ 


(34) 


where  r  is  a  scaling  parameter  governing  the  magnitude  of  a 
point  charge  compared  to  the  one  below  it.  As  charges  of 
only  one  sign  are  used,  y„  in  Eq.  (33)  is  termed  the  “mono¬ 
pole”  model  \^°”°(p,z).  Observe  that  the  first  charge  has 
been  explicitly  pulled  out  of  the  summation  (i.e.,  Xo=  1  and  Zj 
for  7  =  0  is  zero). 

Spindt-type  field  emitters  by  comparison  are  sharper  and 
more  conical.  Therefore,  introduce  a  “dipole”  potential 
defined  by  (no  charge  at  z=0  means  Xo=0) 


n 


Z  +  2  >^j{[p^  +  (z  - 
y=i 


(35) 


which  also  satisfies  y^^P^(0,z„+i)=0,  but  which  produces  a 
sharper  emitter  as  shown  in  Eig.  2  for  r=0.5  and  n  =  6. 

The  trajectories  are  found  by  launching  electrons  with  a 
velocity  hkplm  normal  to  the  surface  defined  by  V^ips^Zs) 
=  0.  The  equations  of  motion  are  numerically  integrated  until 
the  electrons  cross  the  zu  plane,  at  which  point  the  value  of 
P/j  is  found.  Erom  the  relationship  between  p/^  and  kp(ph),  the 
numerical  values  of  the  y’s  are  obtained.  The  value  of  C  is 
obtained  from  the  behavior  of  the  current  density  distribution 
in  the  h  plane. 
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FIG.  3.  (Color  online)  Field  enhancement  factors  as  a  function  of  n  for  the 
mono  and  dipole  emitter  types  for  values  of  r=0.5  and  1.0.  For  large  n,  the 
dipole  held  is  approximately  30%  larger,  even  though  the  apex  radii  are 
nearly  the  same,  for  r=0.5. 


FIG.  5.  (Color  online)  Variation  of  held  along  the  surface  of  a  dipole  ge¬ 
ometry  as  a  function  of  axial  coordinate  p,  showing  the  nature  of  the  held 
drop  off  compared  to  its  apex  value  and  the  quality  of  the  approximation 
cos  6  used  to  derive  Cp^^- 


Field  enhancement,  apex  radii  factors,  and  field 
drop-off 

The  field  enhancement  factor  is  defined  as  the  ratio  of 
the  apex  field  r’tip  to  whereas  the  apex  radius  is  a  ratio  of 
the  first  and  second  derivatives  of  on  axis,  or 


a„{r)  =  - 


dpVnip^z) 

^pV„ip,z) 


p=0 


(36) 


and  they  are  dimensionless.  From  Eqs.  (33)  and  (35),  the 
evaluations  are  straightforward  but  lengthy.  A  comparison  of 
for  r=0.5  and  1.0  for  the  monopole  and  dipole  geometries 
is  shown  in  Fig.  3  and  the  same  for  in  Fig.  4.  Likewise, 
the  quality  of  the  approximation  F{p)=F^^^cos  9{p),  where 
p=a  sin  0,  is  shown  in  Fig.  5,  where  the  line  marked  “Nu¬ 
merical”  is  calculated  using  the  PCM.  The  evaluation  plane 
Zh  is  chosen  such  that  trajectories  from  adjacent  emitters  are 
not  present  in  the  unit  cell,  and  therefore,  Zh-Zo  is  larger  than 
several  tip  radii  but  smaller  than  1. 


model  for  the  “notional”  emission  area^^  [the  ratio  of  total 
current  /^jp  to  current  density  on  axis  /^(F’pp)]  was  developed 
separately but  its  leading  order  behavior  shall  be  used  here. 
We  have 

hipiVs)  =  2'7T-a^g(Ftip)4(Ftip), 


g(F)  = 


0.8F 

[B  +  (l-/c)F]' 


(37) 


The  on-axis  field  F^p  is  obtained  from  the  gate  potential  Vg 
by  F,p  =l3gVg,  where  the  approximation  developed  in  Ref.  35 
is  used  and  is  (for  terms  defined  in  Table  IV) 


77 

In 

a„  1  aA 

—  tan 

8b  + 

_  54(3^  \  (3^  /  _ 

.  (38) 


The  apex  radius  is  adjusted  until  a  theoretical  current  of 
0.5  /uiA  occurs  for  the  same  voltage  as  experimentally  re¬ 
ported.  This  one  data  point  determination  is  not  a  fit  over 
many  data  points:  therefore,  the  degree  to  which  Eq.  (37) 


Relation  to  physical  Spindt-type  emitters 
Emitter  geometry  specification 

The  emitters  characterized  by  Schwoebel  et  al^^  and 
summarized  in  Table  IV  will  be  taken  as  the  base-line  geom¬ 
etry.  The  apex  field  Flip  will  be  inferred  from  FV  data.  A 


FIG.  4.  (Color  online)  Same  as  Fig.  3,  but  for  the  apex  radius:  for  large  n, 
the  estimates  converge. 


TABLE  IV.  Baseline  FEA  unit  cell  parameters.  Parameters  based  on  the  FI 
emitter  shown  in  Figure  1  (I-V  characteristics  of  which  are  shown  in  Figure 
2)  of  Ref.  32.  Ap^  and  Fp^  refer  to  the  least-squares  straight  line  fits  to 
ln(/pN/ y)  ~  ln(ApN) -Fp^/ V  (the  conventional  representation). 


Symbol 

Meaning 

Value 

Unit 

a. 

Tip  apex  radius 

12.4 

nm 

Gate  radius 

0.375 

pm 

Pc 

Emitter  cone  angle 

17 

deg 

dtt 

Tip-to-tip  separation 

1.5 

pm 

dhg 

Base-to-gate  separation 

1.2 

pm 

h 

Base  layer  thickness 

0.18 

pm 

h 

Gate  layer  thickness 

0.75 

pm 

cF 

Work  function 

4.55 

eV 

Ks 

Oxide  dielectric  constant 

3.75 

cr 

Metal  conductivity 

18 

1/(0  pm) 

\ 

Saturn  mode  correction  factor 

1.2 

w 

Square  array  width 

0.1 

cm 

ApN 

Exp.  FN  I(V)  parameter 

8.09 

pA/cm^  V^ 

FpN 

Exp.  FN  I(V)  parameter 

1.72 

kV 

(O 

Gate  sinusoidal  frequency 

277/(5 

1/ps 

S 

Rise  time 

1-50 

ps 
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FIG.  6.  (Color  online)  Comparison  of  analytical  tip  current  model  to  the 
experimental  data  of  Schwoebel  et  al  (Ref.  32)  for  a  single  tip,  in  which 
conditioning  was  shown  to  bring  different  emitters  to  a  comparable  perfor¬ 
mance,  interpreted  as  a  change  in  apex  radius.  Labels  are  those  of 
Schwoebel  et  al.;  lines  are  theory  calculated  using  Eq.  (37). 

reproduces  experimental  curves  is  an  indication  of  its  utility. 
The  tip  radii  found  for  the  curves  labeled  II,  12,  FI,  and 
F2  of  Fig.  2  of  Ref.  32  are  4.17,  10.7,  12.4,  and  12.8  nm, 
respectively.  The  predicted  current-voltage  theory  and  the 
data  of  Schwoebel  et  al.  are  shown  in  Fig.  6.  The  agreement 
indicates  that  Eq.  (37)  successfully  determines  the  permis¬ 
sible  tip  radii:  we  choose  one  corresponding  to  a  “condi¬ 
tioned”  emitter. 

Prediction  of  emittance  for  a  Spindt-type  array 

A  molybdenum  ungated  field  emitter  dipole  PCM  with 
r=0.5,  n  =  7,  and  a  base  radius  of  aQ=l  juum  gives  the  apex 
radius  as  as= a'jaQ=  12.0  nm.  A  background  field  of 
241  MV/m  (the  field  within  the  unit  cell  and  therefore 
higher  than  the  anode-cathode  field)  of  a  gated  FEA  is 
reasonable,^^  corresponds  to  0.241  eV/nm,  and  gives 
Ftip=8.11  eV/nm,  /g(Etip)  =  8.48  X  10^  A/cm^,  and 
=  79.5  nm“^  An  evaluation  plane  at  Zh=2.5  entails  that  rj 
=  1.255.  A  numerical  evaluation  of  the  trajectories  and  an 
extraction  of  y  and  C  from  kp{piJ  and  P{ph)  plots,  respec¬ 
tively,  gives  7i  =  12.7,  73  =  26.8,  and  C=81.4  (compare  to 
estimates  based  on  PSM,  which  give  71  =  8. 15  and  CpsM 
=  65.3,  respectively:  although  of  the  right  magnitude  they 
would  predict  an  emittance  60%  of  its  numerical  value). 

Actual  emitter  radii  in  an  array  follow  a  log-normal 
distribution, making  the  estimation  of  the  number  of 
emitters  operating  for  a  given  gate  voltage  rather  involved. 
An  estimate  of  number  of  tips  emitting  is  reasonably  well 
characterized  by  a  function  A  linearly  dependent  on  gate 
voltage.  Spindt-type  emitters^^  have  shown  current  densities 
comparable  to  2000  A/cm^.  Therefore,  the  current  from  an 
array  is  taken  to  be 

4rray(  Vg)  =  27Ta]g(P^V ^){2N „  +  ^{V ^)J gV g) .  (39) 

The  number  of  active  emitters  is  (2V-H  l)^A(V^)  if  (2V^ 
+  is  number  of  emitters  in  the  array.  It  shall  be  assumed 
that  40%  of  the  emitters  are  active  at  a  gate  voltage  of 
136  V:  therefore,  a  gate  voltage  of  200  V  is  required  to  pro¬ 
duce  a  current  density  of  2000  A/cm^  from  a  square  array 
L=0.2  cm  on  a  side  if  the  tip-to-tip  spacing  is  1.5  ptm.  The 
unit  cell  length  I  is  determined  from  L//=(2V^-h  l)^A(y^), 
and  so  the  emittance  in  the  limit  of  large  is 


I  ^ 

^n,rms  array  “  Tl  /TmTTTTT  ’  v^^) 

^  me  V24CA(y^) 

It  is  therefore  predicted  that  for  a  current  density  of 
2000  A/cm^  from  a  square  array  of  0.2  cm  on  a  side,  the 
emittance  is  23  mm  mrad. 


Comparison  to  experimentai  emittance  measurements 

Comparing  the  theoretical  emittance  formula  to  experi¬ 
mental  data  is  hampered  by  the  scarcity  of  experimental  data 
of  beams  generated  from  field  emitters  (more  often,  beams 
from  photofield  emission  from  arrays^^’"^^  and  single  tips  and 
needles'^ are  reported).  Two  topical  cases  are  Jarvis  et 
al.^^  and  Leeman."^^  Modifications  and  complications  to  the 
theory  are  discussed  for  each. 

Jarvis  et  al.  considered  diamond  pyramidal  field  emitters 
and  find  s^  ^^Qnrvis)  ~  0.97  mm  mrad.  Diamond  and  metal 
emitters  are  categorically  different,  but  semiconductor  emit¬ 
ters  do  have  performance  characteristics  similar  to  metal 
ones  even  though  the  underlying  physics  of  emission 
differs therefore,  as  a  first  approximation,  only  those  pa¬ 
rameters  that  depend  on  layout  and  measurement  conditions 
will  be  altered  from  the  previous  analysis.  Converting  the 
layout  of  Jarvis  et  al.  of  a  tip-to-tip  (pitch)  separation  of 
28  pum  and  the  array  in  the  form  of  a  3  X  24  rectangle  into  an 
equivalent  square,  L  becomes  0.0238  cm.  A  maximum  ex¬ 
traction  field  of  F^=ll  eVZ/xm  was  reported.  Lastly,  the 
emitters  were  conditioned  such  that  A  was  suggested  to  be 
near  unity.  Therefore,  theory  predicts 

I  h  F^maQ 

S«,rms  (Jarvis)  theory  =  Ti — L  “W  — —  =  0.88  mm  mrad, 
me  y  12C 

(41) 

where  the  base  radius  of  the  diamond  emitters  was  taken  to 
be  ^0=2.5  /xm.  As  N  is  not  large  (N=4  corresponds  to  a 
square  array  of  81  emitters,  slightly  larger  than  the  72  of 
Jarvis  et  al),  the  1-h(1/A)  factor  of  Rhee  for  finite  arrays 
suggests  the  theory  value  be  raised  to  1 . 1  mm  mrad. 

Leeman  et  al.  considered  50  000  gated  metal  emitters 
(A^=lll)  emitted  from  a  disk  of  0.1  cm  in  diameter  (L 
=  0.0886  cm)  and  reported  an  emittance  of  2.5  mm  mrad. 
Emitters  that  have  not  been  extensively  conditioned  show  a 
significantly  smaller  A  compared  to  those  that  have,^^  so  that 
on  the  order  of  5%  of  the  emitters  operating  is  expected.  The 
FV  data  of  Leeman  et  al.  in  their  Eig.  2  is  well  modeled  by 
retaining  most  of  the  Spindt-like  parameters  but  changing 
n  =  6  in  the  PCM,  to  give  <2^  =  22  nm  and  E^= 0.215  eV/nm. 
A  comparison  of  Leeman ’s  Eig.  2  data  with  theory  is  shown 
in  Eig.  7.  The  trajectory  simulations  then  give  7i  =  12.1,  73 
=21.4,  and  C=88.2.  An  array  current  of  1.67  mA  at  a  volt¬ 
age  of  195  V  corresponds  to  A(195  V) =0.0412  and  N=23. 
Therefore,  we  shall  project  to  an  array  delivering  1  A/cm^ 
from  an  area  of  0.7854  mm^  at  a  voltage  of  205  V.  The 
smallness  of  A  creates  an  interpretation  problem:  the  active 
emitters  can  be  uniformly  distributed  over  the  array  or,  as  is 
more  likely,  a  hot  spot  of  active  emitters  can  occur.  The  later 
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Gate  Voltage  [V] 

FIG.  7.  (Color  online)  Comparison  of  analytical  array  current  model  to  the 
experimental  data  of  Leeman  et  al.  for  an  array  in  which  on  the  order  of  5% 
of  the  emitters  are  active. 


case  shall  be  modeled  by  assuming  all  the  active  emitters  are 
adjacent  in  a  smaller  square  area.  In  the  “uniform”  case,  I 
=  19  /xm=4.8J^^,  and 

^n,rms  I  spread  out  “  Tl 

^  me  y  12C 

~  6.76  mm  mrad.  (42) 

Conversely,  in  the  “hot  spot”  case,  l=dff  and 


^n,rms  (Lceman)  —  yj  dj^{2N-\- 1) 
me 

—  1.41  mm  mrad.  (43) 

The  measurement  reported  by  Leeman  et  al.  falls  between 
these  limits  and  suggests  that  the  array  is  closer  to  the  hot 
spot  configuration  than  a  uniformly  distributed  one. 


jFomaQ 

\  12C 


3D  comparisons 


considered,  followed  by  performance  values  that  are  more 
indicative  of  how  the  cathodes  are  intended  to  be  used. 


Thermionic  cathodes 

A  circular  dispenser  cathode  of  radius  Pc= 0.5642  cm 
(area=l  cm^)  with  0  =  2.0  eV  subject  to  a  background  field 
of  1  MV/m  and  at  a  temperature  of  r=1200  K  produces  a 
current  of  4=1  A  with  s„  j.jns=L269  mm  mrad.  Larger  total 
currents  are  enabled  by  higher  temperatures  or  larger  cathode 
diameters,  and  greater  charge  in  the  bunch  is  obtained  by 
longer  pulse  durations  and  later  compression  of  the  beam.  If 
the  temperature  is  T(x)  =  (l-\-x)1200  K,  then 

/(x)  =  (l  +  xf 

^n,rms('^)  Vl  ■^^n,rms(^)  •  (44) 

Thus,  to  obtain  100  A  of  current,  the  cathode  radius  can  be 
increased  to  5.64  cm  or  the  temperature  increased  to  1533  K 
(i.e.,  v=0.277),  causing  the  emittance  to  become  12.7  or 
1 .434  mm  mrad,  respectively.  Thermionic  cathodes  cannot 
be  modulated  rapidly  enough  to  produce  pulse  lengths  of 
10-50  ps. 


Photocathode 

A  copper  photocathode  subject  to  a  background  field  of 
50  MV/m,  and  illuminated  with  a  255  nm  laser  of  intensity 
of  10"^  W/cm^  over  a  circular  area  with  radius  of  0.794  cm 
exhibits  a  theoretical  quantum  efficiency  of  0.0246%  and 
produces  a  current  of  4=1  A  with  Sf^  ^^=l.93  mm  mrad.  If 
Pc(v) =0.794(1 -Hv)  cm,  then 

/(x)  =  (l  +x)^4, 


Cathodes  are  operated  under  different  conditions  to  ac¬ 
commodate  their  limitations.  The  beam  drawn  from  the  cath¬ 
ode  experiences  different  acceleration  gradients  depending 
on  the  gun:  for  example,  dc  photoinjector  has  a  smaller  gra¬ 
dient  at  the  photocathode  than  does  a  rf  photoinjector.  Con¬ 
straints  on  the  bunches  themselves  likewise  depend  on  the 
technology.  Therefore,  comparing  the  emittances  of  the  ther¬ 
mal,  photo,  and  field  emission  cathodes  directly  for  same 
cathode  dimensions  and  operating  conditions  would  be  mis¬ 
leading.  Rather,  conditions  reflecting  the  operational  envi¬ 
ronment  typical  of  each  cathode  at  a  performance  level  rel¬ 
evant  to  that  cathode  will  be  considered.  Note  that  emittance 
is  not  the  only  metric.  The  amount  of  charge  in  an  electron 
bunch  is  important  for  accelerators  and  relates  to  current  and 
pulse  duration.  Brightness  goes  as  the  ratio  of  current  to 
emittance  squared  with  high  brightness  be¬ 

ing  desirable.  In  general,  thermionic  cathodes  are  limited  by 
the  source  temperature  (to  avoid  lifetime  issues),  photocath¬ 
odes  are  limited  by  extraction  field  (to  avoid  space  charge 
issues),  and  field  emitters  are  limited  by  maximum  gate  volt¬ 
ages  (to  avoid  breakdown)  or  emission  areas  (to  enable  fast 
switching  times).  Comparisons  therefore  involve  parameters 
which  are  quite  different.  Therefore,  a  baseline  of  1  A  is 


^n,rms('^)  “  ■^)^n,rms(^)  •  (45) 

Modifying  fields  and  intensities  to  match  the  Linac  Coherent 
Light  Source  values  reported  by  Schmerge  et  al  in  which 
a  1  nC  charge  bunch  is  generated  by  a  255  nm  laser  of  radius 
of  1.2  mm  for  approximately  10  ps  subject  to  a  background 
field  of  120  MV/m,  the  current  and  current  density  are  then 
100  A  and  2210  A/cm^,  respectively,  and  the  theoretical 
emittance  is  0.427  mm  mrad.  Photocathodes  can  easily  pro¬ 
duce  pulse  lengths  of  10-50  ps. 


Field  emitter  array 


For  molybdenum  Spindt-type  emitters  with  the  number 
of  emitters  given  by  2A/-I- 1  =  7rp^/ J^^=667  (i.e.,  p^=564  jam) 
with  apex  radii  of  12  nm  and  Table  IV  parameters  otherwise, 
subject  to  a  gate  voltage  of  161.3  V,  then  the  array  current 
density  with  40%  emitting  is  100  A/cm^,  and  the  current  is 
4=1  A  with  mm  mrad.  If  V{x)  =  161.3(l+x)  V, 

then 


lix)  =  1.094  X  10' 


,0+^1 

44.67 -Hx 


exp 


10.11 
1  -\-x 
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^n,rms(-^)  ~  Vl  “I"  ■^^n,rms(^)  •  (46) 

The  maximum  current  density  technologically  practical  from 
gated  field  emitter  arrays  is  likely  on  the  order  of 
2000  A/cm^:  therefore,  if  100  A  of  current  is  required,  then 
the  array  radius  would  be  0.126  cm  and  the  voltage  is  206  V 
(corresponding  to  x= 0.279  and  a  per- tip  current  of  106  /xA 
for  the  active  tips).  The  value  of  e„  rms;(Q)  is  crudely  obtained 
by  scaling  using  e’ ! e  =  {V ! L) Vv/  V  :  taking  L'  =  V^Pc,  V' 
=  y(0),  and  s  =  23  mmmrad,  then  s„  j^s(0)  =  12.8  mmmrad. 
The  theoretical  emittance  of  a  FEA  is  larger  than  for  thermi¬ 
onic  or  photocathodes.  The  demands  per  emitter  may  in¬ 
crease  if  pulse  lengths  of  10-50  ps  from  a  FEA  are  required. 


Issues  of  gated  geometries  and  space  charge 


Useful  field  emission  from  metals  requires  fields  in  ex¬ 
cess  of  5  GV/m  at  the  emission  site,  whereas  fields  in  a  rf 
photoinjector  are  generally  on  the  order  of  10-100  MV/m: 
freestanding  tips  would  have  to  have  field  enhancement  fac¬ 
tors  on  the  order  of  100-1000.  Such  enhancement  would  be 
exceptional  for  metallic  field  emitters:  an  ellipsoidal  height 
to  base  radius  ratio  R  of  16-62  would  be  required  based  on 
an  ellipsoidal  enhancement  factor  given  by  [Eq.  (92a)  of  Ref. 
53] 


2R^-3 

ellipsoid  2[ln(2/?)  -  1] 


(47) 


Even  if  possible,  obtaining  uniform  emission  from  wirelike 
structures  is  difficult  because  of  variation  in  field  enhance¬ 
ment  factors.^"^  To  improve  uniformity,  emitters  are  “condi¬ 
tioned”  in  which  emitters  with  the  largest  field  enhancements 
are  degraded  until  their  values  are  more  representative  of  the 
group. Applications  requiring  not  only  high  current  but 
also  fast  switching  capability  have  therefore  used  gated 
EEAs.^^’^^’^^  The  gate  has  advantageous  and  potentially  prob¬ 
lematic  aspects:  while  the  emitted  beam  is  focused,^^’^^  the 
high  current  density  implies  space  charge  effects,  and  beam 
gating  implies  switching  time  issues  related  to  the  EE  A  cir¬ 
cuit’s  RLC  characteristics. A  brief  and  qualitative  discus¬ 
sion  of  the  impact  of  each  is  given. 


The  gate  and  divergence  of  the  beam 


The  magnitude  of  gate  focusing  can  be  estimated  by 
considering  the  Saturn  model^^’^^  of  a  gated  EE  A,  in  which 
the  emitter  tip  is  replaced  by  a  sphere  holding  a  charge 
the  gate  by  a  coplanar  charged  ring  of  radius  Ug  and  of  equal 
and  opposite  total  charge  -g^,  and  the  gate-anode  field  by  a 
constant  background  field  F^.  The  Saturn  potential  energy  is 
given  by  (in  spherical  coordinates  centered  on  the  tip) 


^)  =  -  For  COS 


61  + 


qQs 

477-So 


1 

r 


CO 

--EC- IV 


r 


>7=0 


(27)! 


r< 

r> 


y 

P2j(cOS  9) 


(48) 


where  Pj{x)  is  a  Legendre  polynomial,  and  where  r<  is  the 


smaller  of  r  or  ag,  and  r>  the  larger.  To  leading  order  in 
(r</r>)^,  the  force  in  the  direction  p  acting  on  an  emitted 
electron  is  given  by 


P  ‘  ^  ^Satum\r<a 


qQs 


sin  9, 


^  .  3a^qQs  ,  , 

P  •  V  7>Satuni  r>a  =  “  - 4  (^9) 

47780^ 

where  r^=z^+p^  (all  terms  of  which  are  dimensioned)  and 
tan  9=plz.  The  presence  of  the  second  term  in  parentheses 
(the  “gate”  term)  is  therefore  seen  to  give  rise  to  a  force 
opposing  an  axial  expansion  of  the  beam  both  within  the  unit 
cell  {r<a)  and  away  from  it  {r>a).  That  is,  the  beam  is 
more  collimated.  However,  the  axial  force  component  will 
both  reduce  the  magnitude  of  the  y’s  and  C,  thereby  mitigat- 
ing  the  impact,  in  line  with  the  arguments  of  O’Shea  that  an 
electrostatic  focusing  structure  will  not  reduce  the 
emittance — or,  alternately,  the  beam  is  already  at  a  waist 
within  the  unit  cell. 


Impact  of  space  charge:  An  estimate 

Eor  high  current  operation,  space  charge  effects  will  ex¬ 
pand  the  beam  and  the  question  is  how  the  magnitude  of 
space  charge  forces  within  the  unit  cell  compares  with  fields 
therein.  Eor  a  crude  estimate,  assume  that  the  apex  field  Fnp 
is  generated  solely  by  g^.  Using  the  parameter  characteristic 
of  the  field  emitter  corresponding  to  line  El  in  Eig.  6  for  a 
current  of  10  julA,  Q^=47TSQa^^^F^^^/ q  —  469q.  Compare  this 
to  the  charge  Ag  existing  within  the  hemisphere  over  the 
tip-gate  region,  approximated  by 

Assuming  /tip=10  /xA  and  Table  IV  parameters,  Ag=16.2^. 
Higher  current  per  tip  will  commensurately  increase  Ag  but 
it  is  clear  that  unless  the  individual  emitters  are  driven  on  the 
order  of  0.3  mA  or  harder,  g^  will  tend  to  dominate  Ag,  that 
is,  the  influence  of  the  gate  exceeds  the  influence  of  the  space 
charge  within  the  unit  cell.  Therefore,  the  emittance  formula 
of  Eq.  (32)  is  approximately  correct. 

The  gate  and  switching  time 

A  thermionic  source  is  run  space-charge  limited  for 
which  the  current  varies  as  a  power  (3/2)  of  the  grid  voltage, 
so  that  the  grid  voltage  must  be  reduced  almost  to  vanishing 
to  turn  off  the  beam.  In  contrast,  field  emission  current  varies 
exponentially  with  gate  voltage  so  that  reducing  the  gate 
voltage  by  a  fraction  of  its  peak  value  suffices. Consider 
the  consequences  of  generating  a  nanocoulomb  charge  bunch 
in  50  ps  from  a  source  required  to  turn  off  completely  in 
10  ps.  The  questions  are:  what  constitutes  “complete  turn¬ 
off’  and  what  consequences  arise  by  insisting  on  a  switching 
time  of  10  ps? 

Eirst,  let  a  be  the  ratio  between  the  minimum  and  maxi¬ 
mum  currents.  Then  the  minimum  voltage  is  the  n>l  limit 
of  the  recursion  relation. 
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FIG.  8.  (Color  online)  RLC  circuit  used  to  model  a  FEA  for  the  evaluation 
of  the  array  length  maximum. 


V. 


BV„ 


n+\ 


B  -H  Vjjiax 


1/  K 


(51) 


With  Vq  ~  Vniax’  convergence  is  achieved  to  better  than  1  % 
after  six  iterations.  We  find  Vjnin/l^max=71-8^  for  C7=0.01 
and  49.9%  for  a=  10“^,  an  advantage  in  reducing  the  drive 
power  required  to  modulate  a  gate.^^  In  short,  the  “extinction 
ratio”  is  excellent. 

Second,  a  FEA  may  be  represented  as  a  bunch  of 
coupled  RLC  units. The  analysis  follows  that  of 
Murphy  and  Kodis^^  to  analyze  the  circuit  suggested  by  Fig. 
8  and  therefore  uses  their  notation.  The  inductance  L,  the 
resistance  R,  and  the  capacitance  C  are  per  unit  length,  and  z 
is  the  direction  in  which  the  signal  causing  the  gate  poten¬ 
tials  to  shift  is  propagating.  Let  the  gate  current  and  potential 
be  denoted  by  Ig  and  Vg,  respectively.  Further,  let  the  sinu¬ 
soidal  frequency  of  the  gate  be  o)  with  the  consequence  that 
the  V^cos(wf)]  curve  can  be  made  to  resemble  a  se¬ 

quence  of  well-spaced  pulses  by  adjusting  the  ratio  of 
.S/Vinax  to  a  number  of  order  of  10.  Using 
=Re[7^(z)^^"T  and  analogously  for  V,  the  transmission  line 
equations 


d^Ig{z)  =  -i(oCVg(z), 


d,Vg{z)  =  -{R  +  io)L)Ig{z)  (52) 

result  and  have  solutions  Ig{z)  =  {-iVo/Z^)sin{l3z)  and  Vg{z) 
=  V^cos{l3z),  where 


P  =  (o^JlC 


R 

1  -  i — 
(oL 


1/2 


(53) 


FIG.  9.  (Color  online)  Variation  of  the  ratio  (l/w),  as  per  Eq.  (60);  as  a 
function  of  array  scaling  factor  5  for  three  different  switching  speeds  (50, 
10,  and  1  ps),  showing  that  the  slower  the  switching  speed,  the  greater  the 
length  of  the  array  can  be. 


scale  factor  such  that  s=l  corresponds  to  dimensions  the  in 
Table  II:  s=2  would  be  twice  as  large,  etc.  The  gate  radius, 
tip-to-tip  separation,  base-to-gate  separation,  and  gate  layer 
thickness  all  scale  linearly  with  i-.  The  apex  radius,  work 
function,  and  cone  angle  in  contrast  are  held  fixed  regardless 
of  scale  size.  As  per  Ref.  35  the  following  scalings  follow: 

/3  j  U  0.080  65  — - ^ - -  -  0.093  47  [ , 

nnm/  [ln[55.4s(2.844  +  s)]  J 

(54) 


5.544 

5’ 


(55) 


;nH\ 

L  1  —  U24.19i, 
^cm/ 


5401 


1.462 


C[8o]  ^  p‘^ln[55.4^(2.844  +  ^)]j’ 


(56) 

(57) 


where  nH=  10  ^  H,  and  units  are  in  the  square  brackets. 
Equation  (57)  is  derived  from 


I-ttI  ^ 


+  (Ug^  +  Ug^)  2  ? 


(58) 


assuming  that  ring  and  cone  capacitances  scale  as  [see  Eq. 
(11b)  of  Ref.  35] 


Q[eo]  = 


3125 


Cg^LsoMifi]  =  10.12s /Bgis), 


Equation  (52)  demonstrates  that  the  gate  voltage  for  each 
emitter  row  changes  with  z,  and  the  emitted  current  changes 
along  the  rows  away  from  where  the  input  signal  is  fed  into 
the  array.  If  an  array  of  width  w  and  length  I  is  to  emit 
uniformly,  then  Evaluating  the  I  that  satisfies  the 

inequality  requires  estimates  of  R,  L,  and  C  for  the  unit  cell. 

The  RLC  unit  cell  characterization  of  a  FEA  of  Calame 

CQ 

et  al.  benefits  from  a  sequence  of  approximations  that  al¬ 
low  for  scaling  arguments  to  be  employed.^^  Let  be  the 


Cgg[soMiii]  =  7.03  U. 

(59) 

Defining  p  =  \pi\,  it  follows 

1  p 

(60) 

For  =  0.1,  the  dependence  of  Eq.  (60)  is  shown  in  Fig.  9  for 
several  values  of  frequency  related  to  the  inverse  rise  time, 
taken  to  be  50,  10,  and  1  ps,  respectively. 


Author  complimentary  copy.  Redistribution  subject  to  AlP  license  or  copyright,  see  http://jap.aip.org/jap/copyright.jsp 


014903-13  Jensen  ef  a/. 


J.  Appl.  Phys.  107,  014903  (2010) 


FIG.  10.  (Color  online)  Ratio  of  array  current  to  single  tip  current  calculated 
by  Eq.  (61)  as  a  function  of  array  scaling  factor  5.  Smaller  emitters  (small  5) 
need  not  be  driven  as  hard  individually  (large  ratio)  for  a  given  total  current. 

The  value  of  I  relates  to  how  hard  the  tips  must  be  driven 
for  a  desired  current.  Letting  the  tip  current  be  given  by  Eq. 
(37),  the  array  current  is 

(61) 

Given  that  wld^^  scales  as  666.1 /s,  then  the  ratio  of  array 
current  to  tip  current  is  shown  in  Fig.  10.  For  the  10  ps  line, 
if  each  tip  is  driven  no  harder  than  0.1  mA/tip  and  the  array 
current  is  (1  nC)/(50  ps)=20  A,  then  /array /^ip  would  have 
to  be  at  least  2 X  10^,  which  is  not  met  (for  the  value  s=l, 
the  ratio  is  a  factor  of  10  too  small).  Consequently,  the  pack¬ 
ing  density,  work  function,  array  width  w,  or  rise/fall  time, 
and  possibly  all,  would  have  to  increase/decrease  according 
to  their  impact. 

CONCLUSION 

Formulas  for  the  emittance  of  a  field  emission  structure 
similar  in  simplicity  to  those  of  thermal  and  photoemission 
have  been  absent  due  to  the  complications  introduced  by 
multidimensionality.  This  work  reports  such  formulas  for 
flat,  single  tip,  and  array  configurations.  First,  a  planar  field 
emission  emittance  equation  was  developed  for  ID  geom¬ 
etries.  Second,  the  consequences  of  an  array  of  emitters  was 
developed.  Third,  a  PCM  was  developed  to  obtain  the  initial 
coordinates  of  emitted  electrons  and  their  relative  proportion 
as  a  function  of  location  on  the  emitter  apex.  Fourth,  the 
method  by  which  the  trajectories  gave  parameters  needed  in 
the  developed  formula  for  emittance  from  an  array  was  pre¬ 
sented.  When  applied  to  Spindt-type  emitters,  the  model  pre¬ 
dicts  that  the  emittance  of  a  square  array  of  0.2  cm  on  a  side 
of  ungated  molybdenum  emitters  to  be  23  mm  mrad. 

Other  comparisons  were  made.  First,  thermal,  photo,  and 
field  emission  sources  were  compared.  Second,  the  magni¬ 
tude  of  space  charge  expansion  forces  within  the  unit  cell 
was  shown  to  be  less  than  the  influence  of  the  gate  and  shall 
be  taken  up  in  a  separate  study.  Third,  the  requirement  that 
the  FEA  be  switched  on  and  off  rapidly  was  related  to  the 
FEA  RLC  parameters,  and  constraints  on  the  field  emitter 
length  were  found,  showing  that  switching  the  arrays  in  un¬ 
der  10  ps  for  performance  conditions  of  interest  herein  may 
be  challenging. 
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